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Supercritical fluids play a significant role in elucidating fundamental aspects of liquid matter 
under extreme conditions. They have been extensively studied at pressures and temperatures 
relevant to various industrial applications. However, much less is known about the structural be¬ 
haviour of supercritical fluids and no structural crossovers have been observed in static compression 
experiments in any temperature and pressure ranges beyond the critical point. The structure of 
supercritical state is currently perceived to be uniform everywhere on the pressure-temperature 
phase diagram, and to change only in a monotonic way even moving around the critical point, not 
only along isotherms or isobars. Conversely, we observe structural crossovers for the first time in a 
deeply supercritical sample through diffraction measurements in a diamond anvil cell and discover 
a new thermodynamic boundary on the pressure-temperature diagram. We explain the existence 
of these crossovers in the framework of the phonon theory of liquids using molecular dynamics 
simulations. The obtained results are of prime importance since it implies a global reconsideration 
of the mere essence of the supercritical phase. Furthermore, this discovery may pave the way to 
new unexpected applications and to the exploration of exotic behaviour of confined fluids relevant 
to geo- and planetary sciences. 
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Statistical mechanics is a very prominent part of physics [I]. In the annals of statistical physics of aggregation 
states, last century and recent decades mark a very vibrant epoch HHH]. A series of successful macroscopic approaches 
suggests that a relatively simple microscopic theory should be capable of yielding realistic phase diagrams. Within 
each phase, the system is uniform in chemical composition and physical state. Critical point occurs under conditions 
of specific values of temperature, pressure and composition, where no phase boundaries exist [9]. As the substance 
approaches critical temperature, the properties of its gas and liquid phases converge, resulting in only one phase 
at and beyond the critical point - a homogeneous supercritical fluid [13 HI]. Recently, it has been experimentally 
established that the supercritical fluids are, in fact, dynamically non-homogeneous in the neighborhood of the critical 
point [T^ . 

The understanding of the supercritical state has recently been revised m- It has been suggested to divide the 
supercritical state beyond the critical point and its neighborhood into two distinct domains by introducing the Frenkel 
line in the framework of the unified phonon-based approach |13j . It has been shown that moving from one domain to 
another in the supercritical state is accompanied by changes in particle dynamics M- The Frenkel line separates these 
two regions beyond the critical point based on changes in phonon excitations, which defines the dynamic crossover 
[Hj. Recently, structural [HKiii and thermodynamic crossovers m associated with the crossing of the Frenkel line 
have been theoretically predicted. Importantly, the predicted structural and thermodynamic crossovers are closely 
related to fundamental changes in phonon states HU US], providing new unexpected connections between elementary 
collective excitations and the structure [HKIS], thermodynamics and scaling laws of supercritical state m- 

In this work, we report results from a diffraction experiment on supercritical argon in a diamond anvil cell (DAC) 
which enables the observation of structural transformations upon crossing the Frenkel line. These transitions corre¬ 
spond to a new thermodynamic boundary on argon pressure-temperature diagram (see Fig.[^. We explain the origin 
of the structural crossovers across the Frenkel line in the framework of the phonon theory of liquids using molecular 
dynamics (MD) simulations, and relate its origin to the change of thermodynamic properties of the supercritical state. 
The experimental results presented here are of great interest to the ongoing effort in elucidating various properties of 
disordered matter pPHSU] - 
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FIG. 1: Argon pressure-temperature phase diagram with the new thermodynamic boundary. The phase diagram 
shows all major states of matter: solid, compressed liquid, supercritical and vapor phases. The inset displays the Argon 
melting line BO], the pressure-temperature tour de experimental path and conditions at which the Frenkel line was observed. 
The behavior of the Frenkel line below 0.4 GPa and 325 K is denoted by the dashed line and is to guide the eye only. The 
position of the Frenkel line at the vicinity (shown by the circle) of the critical point cannot be reliably determined as the critical 
point strongly affects the behavior of all major thermodynamic quantities. At higher P-T conditions, the Frenkel line is nearly 
parallel to the Ar melting line |41| . 


RESULTS 


Mind the transverse phononic gaps 


Here, we introduce the Hamiltonian H = Hq + Hint^ where Hq defines a free theory with no interactions between 
phonons 


OJq<OJ-£) 

and Hint is the interaction term that leads to spontaneous symmetry breaking 


Hint — ^ ^ 


UJq<U)-D 
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g is a multi-index { 51 , 92 , 93 }, wp is the Debye frequency, and the parameter fi takes values 1 or 0. The parameters 
g,X G K"*" are real non-negative couplings, as introduced in the theory of aggregation states [TB]. H" and are the 
collective canonical coordinates and jQ, | = iQqQ-qY^^- The conhgurations Q° and n“ minimise the energy of the 
system and break the SO{3) symmetry to SO{2). Minima of the potential 
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are found to be 
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wf is the Frenkel frequency and defines the lower bound of the oscillation frequency of the atoms or molecules. It can 
be derived from the viscosity r] and shear modulus Goo of a liquid |19) . Excitations of the phonon field around the 
ground state (5“ can be written as 

+<,a = l,2,3 (5) 


where and are the transverse and the longitudinal modes respectively. For a chosen vacuum 
obtain the effective Hamiltonian 
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0<a;q *’*<cdD 


E 

1 1 

-J^q^-q 

+ E 

%^y-q+^y-q) 




- 


Sf\Qq\ we 


( 6 ) 


where I and t stand for the longitudinal and transverse phonon polarizations, respectively. In the framework of the 
above formalism one may derive energy spectra in reciprocal space, varying the system parameters, which in real 
space can be attributed to different states of aggregation such as solids, liquids and gas [TB]. In particular, this 
Hamiltonian predicts that heat capacity at constant volume per particle cy = ^ (fT)N ~ drops down from 
approximately 3 A:b to about 2k^ [nms], where fee is the Boltzmann constant. 

The effective Hamiltonian has a neat property regarding the low-frequency transverse phonon excitations in the 
rigid/compressed liquid regime (see the last term in Eq. [^, i.e. low-frequency wave-packets (long-wavelength limit) 
cannot be propagated in the rigid liquids due to the existence of the transverse phononic gaps in the spectrum 
(0 ^ wf < w*’* < wd). The inability to support low-frequency transverse elementary collective excitations (low- 
energy cutoff) is a manifestation of the absence of the long-range order structural correlations in liquids (see DAC 
experiments) and vice versa. An increase in temperature leads to the disappearance of both the high-frequency 
transverse phonon modes ( wf —> wd, hence, cy'- 3A:b —> 2/i:B) [HI 112 03] and progressively the medium-range order 
pair correlations [mill]. cy = 2/cb ( wf = wd, see Eq. is the new thermodynamic limit ( dubbed here the Frenkel 
line thermodynamic limit) along with other well-known cy = (the Dulong-Petit law) and cy = |fcB (the ideal gas) 
thermodynamic limits which are also covered by Eq. Therefore, crossing the Frenkel line results in fundamental 
changes of pair structural correlations both in reciprocal and real spaces [m [Ti|- as well as in the thermodynamics 
(heat capacity at constant volume cy), scaling laws |17j . and phonon states in the supercritical matter jl81119j . 

Alternatively, the liquid energy can be calculated on the basis of pair correlations 


E = -Nk^T + 271N p 


r^u(r)g(r) 


( 7 ) 


where p is the density, u(r) is the interatomic potential, and g(r) is the pair distribution function which can also be 
represented as a Fourier transform of the static structure factor S(q) [15] . 

The first approach (see Eq. takes into account the phonon contributions into the internal energy of a liquid. As 

temperature increases (wf —> wd) the internal energy and the heat capacity approach the Frenkel line thermodynamic 
limit meaning that cy = ^ 2fcB. On the other hand, from the second equivalent approach (see Eq. 7), the 

internal energy with temperature increase loses the contribution from the medium range pair correlations whic 


n can 


be evidenced from disappearance of the 2”*^ g{r) or S{q) peaks. Thus, cy = ^ 2A:b. Therefore, dynamic, 

structural and thermodynamic crossovers are interconnected within the liquid and the supercritical state and governed 
by temperature variations, hence, negligibly influenced by the pressure in terms of the described analysis. 

It is noteworthy, that any approach what bears a ’’small parameter” and makes it possible to expand the energy 
of the system into the series, for example such as virial expansion, Meyer expansion, Percus-Yevick approach, hard- 
sphere model and other perturbation approaches are capable to describe weakly interacting systems, by the definition. 
As a result, these approaches fail in describing the systems with strong interactions like compressed liquids. Hence, 
these methods are unable to predict and describe the structural crossover between compressed and non-compressed 
liquid regimes which is described above and the subject of this work. 


Diamond anvil cell experiments 

The combination of synchrotron X-ray facilities with high-pressure methods provides new experimental tools for 
probing the structure and dynamics of materials. In our experiment, we performed high pressure/high temperature 
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FIG. 2: Experimental reciprocal space evolution of pair correlations across the Frenkel line.(A) The measured 
static structure factor S{q) on X-ray diffraction experiments at different pressure-temperature conditions. (B, C) First and 
second diffraction S{q) peaks pressure-temperature position variation. Labels indicate the pressure in GPa units. Insets show 
peaks pressure-temperature heights variation illustrating non-uniform behaviour, hitherto unanticipated |10| . (D) Comparison 
of the experimentally obtained S{q) with the S{q) calculated from MD simulations. The inset shows the g{r) derived from the 
MD simulations. 


X-ray diffraction (XRD) measurements using a BX90 (DAC) [Hj at GSECARS beam line 13-ID-D of the Advanced 
Photon Source (APS), Argonne National Laboratory. The DAC was used in combination with tungsten-carbide seats 
and full diamond anvils with a 500 /rm culet size. 250 /rm-thick rhenium gasket was pre-indented to a thickness of 
about 40 /im. A hole with a diameter of about 120 /rm was drilled in the middle of the pre-indented area. Conventional 
resistive heating was used to heat the sample up to 500 K. The temperature was measured by two thermocouples 
attached to the diamond surface at the vicinity of the culet. The '*°Ar was loaded using a COMPRES/GSECARS 
gas-loading system at APS [iS] up to initial pressure of 1 GPa. A ruby sphere was used for the pressure calibration 
[46] . Given the Ar loading pressure, the thickness of the gaskets pre-indented area and the diameter of the drilled hole, 
the sample would amount to approximately 10^® Ar atoms. The sample was measured at the photon energy of 37.077 
keV; the XRD patterns were recorded by a MAR-165 CCD camera with 79x79 /rm^ pixel size with the exposure time 
of 10 seconds per XRD pattern. Following every temperature change, the DAC was allowed to equilibrate for, at least, 
5 minutes before the XRD pattern was collected. The raw XRD patterns were further processed using the Dioptas 
software package. The background measured from the empty cell was subtracted from each spectrum. Each S{q) 
curve was normalized to the Ar compressibility limit (see, e.g. page 31, |47| 1 for a given pressure and temperature. 
The 5'(0) values are obtained from independent thermodynamic compressibility data (see the NIST database). 

In the present definitions of supercritical state [M E], even moving around the critical point nni HI] on the 
pressure-temperature diagram (not only along isotherms or isobars) can be considered appropriate to detect possible 
structural non-uniformities based on pair correlations analysis [101 [II]. Selected S{q) curves at different P-T conditions 

















































5 


for wavenumbers q up to 6.5 A“^ are presented in Fig. [^. In Fig. [^b,c), we present the position variation of first 
and second S{q) peaks as a function of temperature. Insets show the S{q) peak height as a function of temperature 
illustrating the non-uniform behaviour consistent with previous theoretically predicted results [mile]. 

We would like to stress that the evolution of the S{q) peak positions (see Figs. ib, c)) as a function of just pressure 
is featurelessly monotonic and has no implications for the observation of the Frenkel line. As we stated above, the 
thermodynamics of supercritical fluids shall be studied as a function of temperature with a little regard to the pressure 
deviations along the experimental path. Indeed, as can be seen from the previous studies [JS] the dependence of the 
first and the second Ar S{q) peak positions as a function of pressure at the fixed temperature (isothermal scan) is 
essentially linear. Clearly, the mere linear variation of the S{q) peak positions as function of pressure cannot affect 
the complex evolution of the peak positions presented in the Figs. Hb, c)(also see Supplementary Materials). 

Therefore, the structural crossover can only be observed as a function of temperature (in contrast to previous works 
[48| where isothermal scans were studied and no crossovers detected) and explained within the framework of the 
phonon theory of liquids where temperature (not the pressure) is the key variable. 


MD simulations at experimental conditions 

In order to analyze the pair correlations in real space, see Fig. [^a,b), we performed MD simulations at experimental 
conditions. In addition, we compared the experimentally determined S{q) with S{q) derived from MD simulations 
showing good agreement, see Fig. and consistent temperature evolution, see Figs. [lb,c)-§c ,d). The structure 
factor S{q) can be defined as 

S{q) = 1-1- Airg / drr^ - — ('?(''’) ~ 1) (8) 

Jo Qr 

where the g{r) is the pair distribution function [S], it describes the distribution of distances between pairs of particles 
contained within a given volume, and Rmax is the distance cutoff parameter set to 20 A which we have found to be 
sufficient to converge the integral. We have used LAMMPS simulation code to run a Lennard-Jones (LJ, e/kB=119.8 
K, a =3.405) fluid htted to Ar properties [JS] with 32678 atoms in the isothermal-isobaric (NPT) ensemble. We have 
used 3000 processors of a high-throughput cluster with a runtime over 300 picoseconds. 

In Fig. [^we observe the alternation of regimes, which is a result of power laws change, both for the first and the 
second peak positions and heights resulting from distinct changes between medium- and short-range order correlations 
[lain]. Different regimes are manifested in different slopes of peak positions and their heights as a function of 
temperature (red and blue dashed lines are guide for the eyes only). The current interchange of the medium- and 
short-range order correlations is reflected in the last term in Eq. The presence or the absence of the last term 
in the Hamiltonian describes structure and thermodynamics above and below the Frenkel line, respectively, which 
can be referred to as different regimes. Remarkably, experimental evolution of pair correlations (see Figj^ also 
illustrates non-uniform behavior as simulated does (see Fig. [^. The non-uniform behaviour of pair correlations both 
in reciprocal and in real spaces is closely related to the differences in the relaxation processes, phonon excitations and 
thermodynamics above and below the Frenkel line [IZIIIH]- 

Although the structural order of a fluid is usually enhanced by isothermal compression or isochoric cooling, a few 
notable systems show the opposite behaviours. Specifically, for example for H 2 O m. increasing the system density 
can disrupt the structure of liquid-like fluids, while lowering temperature or strengthening attractive interactions can 
weaken the correlations of fluids with short-range attractions. Structural order is a particularly useful measure of 
the observed changes crossing the Frenkel line. It also correlates strongly with the static structure factor and the 
pair distribution function behaviour which provides insight into the structural crossover in the supercritical state. At 
elevated temperatures (> 1000 K), decrease of the first (both S{q) and g{r)) peak and the near disappearance of the 
second and third peaks imply the medium-range order correlations are no longer present. Such high temperatures 
are not readily accessible in DAC experiments but can easily be studied using MD simulations. Nevertheless, the 
experimental evidence for this behavior is clearly seen as a decrease of the 1®* and 2"^^ S{q) peak heights above the 
T ^ 410 K as shown in Fig. [2];b,c) (see insets). The peak height decrease is also presented in Fig. [I] and Fig. [| 
The temperature variation of the 1®* and 2”®* S{q) peak positions is also similar in the DAC experiment and MD 
simulations in the temperature range above T ~ 410 K. However, the structural crossover at elevated temperatures 
cannot be observed experimentally in DAC measurements as the transition takes place at ^ 1500 K. 
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FIG. 3: Real space evolution of pair correlations across the Frenkel line as derived from MD simulations.(A, 

B) and (C, D) first and second g{r) and S{q) peaks pressure-temperature position variation derived from the MD simulations 
at experimental conditions, respectively. (A, C) labels indicate the pressure in GPa units. Insets show g{r) and S{q) peaks 
pressure-temperature heights variation exhibiting heterogeneous behaviour, hitherto unanticipated [IS]. The dashed lines are 
to guide the eye only. 
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FIG. 4: MD simulations: structure temperature variation in reciprocal space. (A) Evolution of static structure 
factor S{q) showing the disappearance of the medium-range correlations with increasing temperature. (B,C ) First and second 
S{q) peaks heights temperature evolution. Insets show S{q) peaks and their derivatives plotted in linear scale. Arrows in 
the insets indicate the corresponding Y-axis for each curve. (D) The position of the E* and 2"'* S{q) peaks as a function of 
temperature. Arrows indicate the corresponding Y-axis for each curve. (B, C, D) The graphs highlight two different regimes 
with two unique power laws (denoted by red solid straight lines) governing S{q) peaks positions and heights resulting from the 
crossing the Frenkel line. 


MD simulations at extended temperature range 

To analyze the behaviour of pair correlations in detail above and below the Frenkel line we run the MD simulations 
within very wide temperature range, which is not accessible for DAC experiments. Here, we have simulated a one- 
component Lennard-Jones (LJ, e =0.994, a =3.405) fluid fitted to Ar properties [5T]. The system is composed 
of a constant-volume (NVE) ensemble of 32000 atoms. Simulations were performed over thewide temperature range 
extending well into the supercritical region (see Fig.|^. The temperature range used is located approximately between 
3Tc and 167Tc, where T^. is the critical temperature of Ar, T^. ~I50 K or 1.3 in LJ units. The system was equilibrated 
at constant temperature. The simulated density, 1880 kg/m^ (1.05 in LJ units), corresponds to approximately three 
times the critical density of Ar. A typical MD simulation lasted about 50 picoseconds (ps), where the pair distribution 
function g{r), static structure factor S{q) and self-diffusion were averaged over the last 20 ps of simulation, after 30 
ps of equilibration. A characteristic relaxation time of any liquid, including Ar, is about 0.1-0.3 ps, therefore the MD 
simulations significantly exceeded the mean collision time. Simulations were performed at 100 temperature points 
within the temperature range of interest. The temperature variation of these quantities is shown in Figs.[4][^ 

In the limiting case of no interaction, the system behaves as an ideal gas and the structure factor is completely 
featureless: S{q) — 1. There is no correlation between the positions rj and of different particles (they are 
independent random variables) so the sum over the off-diagonal terms in equation 



( 9 ) 









































FIG. 5: MD simulations: structure temperature variation in real space. (A) Evolution of mean force potential, (B) 
diffusion and (C) its first derivative and first coordination shell at different temperatures. (D) (o) = 0.5 x giMax{r)x I”** peak 
position is plotted in Log-Log scale and highlights two different regimes with two unique power laws (represented by red solid 
straight lines). Inset shows the first coordination shell plot in Log-Log scale. Arrows indicate the corresponding Y-axis for each 
curve. 


vanishes, i.e. 


(exp[-fq(rj - r^)]) = (exp(-'jqrj)) (exp(fqrfe)) = 0 (10) 

Even for interacting particles at high magnitudes of the scattering vector q, the structure factor goes to 1. This result 
follows from equation 


S{q) = 1 + p [ dre (11) 

Jv 

since S{q) — 1 is the Fourier transform of the g{r) and thus goes to zero for high values of the argument q. This 
reasoning does not hold for a perfect crystal, where the distribution function exhibits infinitely peaks. To examine the 
temperature changes of S{q) in more detail (see Fig.|^), we show temperature variation of the heights and positions 
of the first and second peaks of S{q) in Fig. Bb -d). There is a steep decrease of both peaks at high temperature. 
Fig. Qb,c), followed by their flattening at high temperature, with the crossover between the two regimes taking place 
around 1500 K. The crossover is visible in the insets of Fig. BbjC) where we plot the temperature derivative of the 
heights of both peaks. These plots clearly show two structural regimes corresponding to the fast and slow changes of 
the S{q) peaks and their positions. The structural crossover in the reciprocal space of supercritical state at elevated 
temperatures which we observe in Figs. |^b-d), supports the experimental results from the main text. It should 
be noted, that the observed crossover is not related to the melting line, as the S{q) profile of the Ar solid phase 
(well-ordered phase) would exhibit very clear sharp peaks throughout the q-range covered by MD simulations. In 
contrast, the 1®* and the 2"'^ S(q) peaks become wider and weaken with temperature increase signifying the strong 
disorder. 

The structural crossover is further evidenced by the calculation of mean force potential (MFP), self-diffusion and 
first coordination shell parameters that include the first shell coordination number, the height and the hrst g{r) peak 
position at different temperatures (see Fig. [^. The temperature variations of the system properties are manifested 
in both real (Fig. and reciprocal (Fig. |^b-d)) spaces. 
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Various methods have been proposed for calculating the MFP. The simplest representation of the MFP is the 
separation r between two particles as the reaction coordinate. The MFP (see Fig. [^) is related to the g{r) using the 
following expression for the Helmholtz free energy (52] 

w =-kBT\n[g{r)] (12) 


where fee is the Boltzmann constant. We use Einstein’s formula to calculate self-diffusion from the mean square 
distance (MSD) travelled by a certain particle over a certain time interval. At the limit of infinite observation time, 
self-diffusion in terms of MSD [SSI becomes 




(13) 


where r) is displacement vector of the i-th atom at t time and the term [fi{t) — h)(0)]^ is the MSD. 


DISCUSSION 

The relationship between system properties in the real and reciprocal spaces is of prime interest in condensed matter 
physics. It is widely recognised that such a relationship may exist in some classes of systems but not in others. In this 
work, the DAC experiments and MD simulations show that not only the supercritical state is physically non-uniform 
but also that manifestations of the non-uniformity are observed in both real and reciprocal spaces (S'(g) and g{r) peaks 
exhibit a ’’zigzag” behavior). We relate the structural crossovers in the supercritical state to changes in elementary 
collective excitations, power laws and thermodynamics in the crossing of the Frenkel line. Crossing the Frenkel line 
corresponds to a quantitative change of the supercritical fluid’s atomic structure and the transition of the substance 
from the compressed liquid structure to the non-compressed gas-like structure (see Figures [2][^ and discussion below). 

The structural crossover takes place continuously when the liquid relaxation time r (the average time between 
two consecutive atomic jumps at one point in space |I81 II9j l approaches its minimal value td, the Debye vibrational 
period, upon which the system loses its ability to support high-frequency propagating shear modes with to > '^ and 
behaves like a gas (EHni- When all shear modes are lost, only the longitudinal mode remains in the system and the 
heat capacity at constant volume becomes cy = 2fcB per particle [mui]. This result can be easily obtained from 
the phonon theory of liquids in the classical limit m- Therefore, the observed thermodynamic boundary in the DAC 
experiments is also closely related to both dynamic and thermodynamic crossovers existing in the supercritical state 

[lain]. 

The thermodynamic boundary discovered in this work has several possible implications to astrobiology and the 
existence of alien life with the impact of extreme conditions on biomolecules |54j . Astonishingly, a number of species 
of bacteria are tolerant of supercritical CO 2 and can survive under severe conditions of pressure and temperature [55] . 
Recently, we have studied structural properties of the supercritical CO 2 |I6j where we provided the evidence for the 
existence of persistent medium-range order correlations that make supercritical CO 2 nonuniform and heterogeneous 
on an intermediate length scale. In the first shell of the CO 2 cluster both carbon and oxygen atoms experience 
gas-like type correlations with short-range order interactions while within the second shell, oxygen atoms essentially 
exhibit a liquid-like type of correlations due to localization of transverse-like wave-packets. Atoms inside the nearest- 
neighbor heterogeneity shell play a catalytic role providing a mechanism for diffusion on an intermediate length 
scale. Extraterrestrial organisms might use these peculiar structural and thermodynamic advantageous properties 
of the supercritical CO 2 to survive biologically. Other possibilities include exoplanets with 2-5 times the mass of 
the Earth with stronger gravitational pulls (super-Earths) and, thus, having supercritical atmospheres and/or oceans 
m- Therefore, further theoretical and experimental study of thermodynamics, dynamics and structure of various 
supercritical fluids and supercritical carbon dioxide in particular are in a strong demand. We believe that this 
discovery will boost the industrial use of supercritical fluids in more efficient way, assist us to search for other 
alternative conditions for an extraterrestrial fine-tuned life and also can lead to greater understanding in another 
disordered systems such as glasses and granular materials. 
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